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Abstract 

We examine the equilibrium and stability of an elastocapillary sys¬ 
tem to model drying-induced structural failures. The model comprises a 
circular elastic membrane with a hole at the center that is deformed by 
the capillary pressure of simply connected and doubly connected menisci. 
Using variational and spectral methods, stability is related to the slope of 
equilibrium branches in the liquid content versus pressure diagram for the 
constrained and unconstrained problems. The second-variation spectra 
are separately determined for the membrane and meniscus, showing that 
the membrane out-of-plane spectrum and the in-plane spectrum at large 
elatocapillary numbers are both positive, so that only meniscus perturba¬ 
tions can cause instability. At small elastocapillary numbers, the in-plane 
spectrum has a negative eigenvalue, inducing wrinkling instabilities in thin 
membranes. In contrast, the smallest eigenvalue of the meniscus spectrum 
always changes sign at a pressure turning point where stability exchange 
occurs in the unconstrained problem. We also examine configurations 
in which the meniscus and membrane are individually stable, while the 
elastocapillary system as a whole is not; this emphasizes the connection 
between stability and the coupling of elastic and capillary forces. 


1 Introduction 


Elastic deformations induced by capillary forces have been identified as leading 
causes of pattern collapse in miniature electronic devices and sensors [l^, [T3| . 
These microstructures are more prone to collapse when miniaturized because 
adhesion and capillary forces become comparable to elastic forces when there 


is a high surface area to volume ratio [ll|. Capillary driven collapse poses 
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major problems for micro-fabrication techniques that are based on wet etch¬ 
ing. In particular, microelectronic systems are commonly fabricated through 
wet lithography where structures often experience permanent deformation and 
stiction upon drying, significantly limiting the design and operating conditions 


34|. 


ElastocapiUary systems have been extensively studied over the past two 
decades 0)[l6, 3, 25, 3 3 . Aggregat ion, coalescence, and self-assembly of fila¬ 
ments and flexible fibres 13 321, failure of microelectronic devices 2^ 331, 


failure of microelectronic devices 

and capillary wrinkling of elastic membranes (^ . 1^ are applications where 
system configurations and structures are determined by elastic-capillary force 
interactions. The foregoing studies are mostly concerned with systems where 
equilibrium configurations are always stable (or assumed to be). However, me¬ 
chanical stability is central to applications in which preventing structural failure 
upon drying is crucial. 

Recently, a few studies have focused on the mechanical stability of elastocap- 
illary systems. Giomi & Mahadevan [l9| examined the equilibrium and stability 
of minimal surfaces spanning deformable frames. Subjecting a circular frame 
to spatial perturbations, they approximated instability modes and the critical 
elastocapillary numbers corresponding to the primary and secon dary buckling 
of the frame into elliptical and twisted structures. Taroni & Vella [38[ identified 
multiple equilibria in an elastocapillary system related to the aggregation of 
paint-brush bristles where the stable solutions for a given liquid content were 
determined through a t emp oral stability and dynamic analysis. 

Mastrangelo & Hsu [25j took a different approach to determine stability in 
elastocapillary systems. Their approach hinges on a pervasive theory, known as 
catastrophe theory Q in nonlinear dynamics, stating that stability exchanges 
only occur at folds and branch points on equilibrium branches (^ . While this 
has not been generally proved for all mechanical systems, the idea has been 
extensively examined for purely capillary MiEi” 1^ . lial and purely elastic 
problems. In this context, Maddocks Q established a theory for systems 
where equilibria are described by a continuous functional of a single function 
with prescribed boundary conditions (from a Hilbert space). This theory relates 
the stability of constrained and unconstrained variational problems to the shape 
of equilibrium branches with no branch point where stability exchanges occur 
only at simple folds. 

Determining the stability of an elastic structure deformed by the Laplace 
pressure or contact line force of a meniscus is more challenging than determining 
its equilibria. Elastic and capillary parts for equilibrium states can be decoupled 
and determined separately by imposing the proper boundary conditions where 
the meniscus and structure meet. However, the stability of elastic and capillary 
parts alone is not sufficient to deduce the stability of elastocapillary systems for 
which the control-parameter role is particularly important. 

In the present work, we study an elastocapillary problem where confor¬ 
mations are controlled by the liquid content, similarly to Kwon et al. 22|. 


This is analogous to the problem of disconnected free surfaces considered by 
Slobozhanin |36l | . Complexities in elastocapillary problems arise because menis- 
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cus perturbations are neither pressure controlled nor volume controlled, as in 
purely capillary problems [1, [i^. Moreover, elastic structures and menisci in 
many practical elastocapillary systems have free boundaries @,[3 US [Hi! which 
considerably complicate the stability analysis. Vogel [H| highlights two major 
difficulties for analyzing the quadratic forms arising from the second variation 
of systems with free boundaries: (i) The function space of perturbations is not 
necessarily a symmetric Hilbert space . Instead, the quadratic forms are 
naturally expressed in and an additional analysis is required to link the 
arising operators to the corresponding operators in a symmetric spac4i|. 
(ii) Perturbed surfaces resulting from normal variations of an equilibrium sur¬ 
face are not generally guaranteed to satisfy the boundary conditions at the free 
boundaries. 

In this paper, we examine the elastocapillary system shown in Fig. [T] and 
relate stability to the slope of equilibrium branches in pressure versus volume 
diagrams, similarly to Maddocks |24| . This system is a model for drying-induced 
structural failures arising in practical applications, such as the stiction of micro- 
machined sensors and collapse of wood fibres upon drying. It comprises a cir¬ 
cular elastic membrane, with a hole at the center, anchored above a rigid plate, 
trapping a prescribed volume of liquid. We examine membrane deformations 
caused by a meniscus at the hole as the liquid is slowly removed. Our approach 
is variational, so that linear stability is determined by the sign of the second 
variation. We demonstrate that there are configurations in which the menis¬ 
cus and membrane are individually stable, while the elastocapillary system as 
a whole is not. This emphasizes the significance of instabilities arising from the 
coupling of elastic and capillary forces. This result can be interpreted as the 
equivalent of the Weierstrass-Erdmann condition 0 for the second variation, 
and it is relevant to applications where extrema are represented by non-smooth 
functions, such as for elastocapillary systems, threshold phenomena 1^, and 
data visualization 1281 . 


2 Formulation 

We consider an elastocapillary model comprising a circular elastic membrane 
with a hole at the center supported on the sidewall of a cylindrical cavity with 
rigid walls, trapping a liquid volume vi below the membrane and air volume 
Vg between the bounding surface (dashed line in Fig. [TJa)) and membrane, as 
shown in Fig.[TJ The cavity is open to the atmosphere from the top. A meniscus 
forms at the hole as the liquid is removed, resulting in a difference between the 
liquid pressure pi and atmospheric pressure Pg , which causes the membrane to 
deform. Here, the membrane radius R, hole radius i?o, and cylinder hight H are 
the model length scales that control the interplay between elastic and capillary 
forces. To determine the equilibria at a given vi, we consider an imaginary 

^Note that denotes Sobolev spaces equipped with the Euclidean norm [J|. 

Moreover, throughout the paper, ‘symmetric space’ refers to a space in which all bilinear 
forms are symmetric. 
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Figure 1: Elastocapillary model; (a) schematic showing simply connected menis¬ 
cus (top), doubly connected meniscus (bottom), transition from simply to dou¬ 
bly connected meniscus (middle), and (b) contact angles. 


bounding surface that covers the cavity from the top. The system is completely 
isolated from the surrounding by the bounding surface and cylinder walls. The 
meniscus is initially a bubble, which can bridge the gap upon contact with 
the plate at the bottom of the cylinder, forming a free contact line with the 
plate. Assuming that all the dimensions are small compared to the capillary 
length, the gravity force is neglected. The membrane and meniscus are assumed 
axisymmetric in equilibrium and perturbed configurations. 

2.1 Variational principle 

Following the development of Neumann et al. , we apply a variational prin¬ 
ciple to determine stable equilibria for the system depicted in Fig. [T] Noting 
that the liquid volume is the control parameter in drying, the total internal 
energy is to be minimized subject to vi = const., maintaining fixed total en¬ 
tropy and mass. Here, the grand canonical potential is a suitable free-energy 
representation because it restricts the minimization to states that are already in 
thermal (constant uniform temperature) and chemical (constant uniform chem¬ 
ical potential) equilibrium. The grand canonical potentials for bulk phases and 
interfaces, respectively, are 

_ y(a) _ 2 ig(a) _ ^^(a) _ a = gl,sl,sg (2) 

with 7 , w, u, s, and p the surface tension, specific grand canonical potential, 
specific internal energy, specific entropy, and density of the respective phase. 
The superscripts (v) and (a) denote volume density and area density for bulk 
phases and interfaces. Note that the temperature T and chemical potential p 
can be regarded as the Lagrange multipliers associated with the entropy and 
mass in the foregoing constrained minimization of the total internal energy. 
Hence, stable equilibria minimize 

Et = (3) 

where F^- are interfacial surface areas. Here, the membrane strain energy 
Ol”") is separately incorporated into the total energy Et to account for vari¬ 
able and anisotropic stresses. Neglecting the bending energy, we only consider 
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the stretching part of the elastic strain energy in von Karman’s theory for mod¬ 
erately large deflections 

— f (^NrrSrr + NttStt)dA, (4) 

where Fmo, Nu, and eu are the membrane in the referential configuration (un¬ 
deflected state), axial forces, and nonlinear strains [i^. Substituting Eqs. ([T]), 
0, and (III) into Eq. ([3]) and omitting additive and multiplicative constants that 
do not affect the minimization, 

Et[r,u,w,P] = U[r,u,w] — PJ[r,u,w], x x x R —>• M (5) 


is the functional to be minimized subject to vi = const., where 


17[r, u, in] 



F{z, r, r')dz 



u, w, u', w')drp + -^ilsg 


7si), (6) 


ph pR 

J[r, M, ic]= / K{z,r,r')dz+ / M(rp, it, w, m', ic'jdcp-|- 

j Zq J Rqo 


Rlh 


with integrands 


F(z,r, r') = 7gjr\/r + 

p 2 

K{z,r,r') = - — , 


C 


2i'uu' vuw'^ 


G{rp, u, w, It', w') = TT H-1-1" “9 + 


It w 
2 ■ 4 


Ip Ip • p 

M{rp,u,w,u',w') = uH + {jp + u)H -|- {vp + u){H -|- w)u'. 


(7) 

( 8 ) 
(9) 

( 10 ) 

( 11 ) 


Note that C, v, it, w, i?i, and i?oo are the membrane axial rigidity, Poisson 
ratio, membrane in-plane displacemenl[^, membrane deflection, radius of the 
meniscus contact line with the plate, and hole radius in the referential configu¬ 
ration. Here, primes denote derivatives with respect to the function argument, 
P = Pi — Pg can be regarded as the Lagrange multiplier associated with the 
constant vi constraint, and the menisci are represented by r{z). The membrane 
deformations are represented by u{rp) and wirp), where rp is the radial coor¬ 
dinate in the referential configuration. When the meniscus is a bubble (simply 
connected), the last term in Eq. (|6|) is zero and zq = where the meniscus in¬ 
tersects the symmetry axis. When the meniscus is a bridge (doubly connected), 
zo = 0, where the free contact line rests on the plate. 

^Not to be confused with the specific internal energy in Eqs. o and ©■ 


5 







2.2 Equilibrium from first variation 

To construct the increment of Et in Eq. ([S]) with respect to axisymmetric per¬ 
turbations, perturbed states are represented by 


z(z) = z + r 7 i(z)e -f r] 2 {z)e'^, 

■■ [6,(i] -t K, 

(12) 

r(z) = f(z) -1- ^i{z)e -1- 6(2)e^i 

7,6,6 : [zo,h] K, 

(13) 

u{rp) = u{rp) -1- ^i(rp)e -|- ^2(?'p)e^, 

fi, 4>1t4>2 ■ [77oo, R] R, 

(14) 

w{rp) = w{rp) + tpi{rp)e + tj}2{rp)e'^, 

w,'0i,•02 : [7?oo,77] -t R, 

(15) 


accounting for the linear and nonlinear parts of the increment when the meniscus 
is a bridge. Note that the form of the functionals in Eqs. dH]) and (Q demands 
f', u', w' to be continuous, so r,u,w € L^DC^. Here, equilibrium and perturbed 
states are denoted by hatted and unhatted variables, respectively. Equation m 
is particularly important, because it admits perturbations that can displace the 
position of the bubble apex and hole edge along the z-axis. We impose a si mply 
supported boundary condition for the membrane at rp = i? where u,w = 0 |4fl| . 
resulting in 

MR) = MR) = MR) = MR) = 0- (16) 

The meniscus is assumed to be pinned to the hole edge at = i?oo, where 
= ?'(-Roo)|r,„ and z{h)\r^i = z(i?oo)|r„, furnishing 

Rq = Roo + u{Rqq) , = (^i(i?oo), ^2{h) = 4'2 {Roq), (17) 

h = H + w{Roo), 77i(/i) ='0i(i?oo), V 2 {h) ='(p 2 {Roo)- (18) 

Moreover, the meridian curve intersects the symmetry axis at z = i, where it 
can only move vertically when the meniscus is a bubble, so 

Ci(f) = 6(f) = 0, ??i(f),? 72 (f) = finite, (19) 

whereas the contact line with the plate can only move horizontally at z = 0 
when the meniscus is a bridge, so 

6(0),6(0) =finite, r;i(0)= 772(0) = 0. (20) 

Since the domain of r(z) is variable in Eqs. m and 0, the functional 
variations are properly rep resented with respect to the barred component of ^ 
(see Gelfand & Fomin [l^ for details) 

6=6+?7i7', (21) 

^2 = ^2 + V2r + . ( 22 ) 


6 


Substituting Eqs. (jl2p - (jl5p into Eq. the first variation of Et with respect 
to an equilibrium state is 


5Et 

e 


+ 


Pr'\h — P Roh — Gu'Irqo + 


- PJ[w),^l) 

Uh) 




h ^ Gu 


Roo]VlCh)+ Rlhsg - Isl) - Er-'lo Cl(0), (23) 


where all the functional integrands are evaluated at the equilibrium. Here, 
primes operating on functionals denote the first Frechet derivative Q with re¬ 
spect to the function in the subscript, and (•,•) is the inner product over the 
domain of the respective function. The last term in Eq. (1231) must be replaced 
with [f'Fr' — F]^7 ]i{£) when the meniscus is a bubble. Equilibria are the sta¬ 
tionary points of the total energy where SEt = 0 for arbitrary ^i, ^i, and ' 01 , 
requiring 

Ulf) - PJ[f) = ^ ( 24 ) 

U[u) - PJ[u) =Gu- PMu - {Gu' - PM^,) = 0, (25) 

U[^) - PJ[^) =G^- PM^ - ^ {G^, - PM^,) = 0 , (26) 

with each boundary term in square brackets equal to zero. Substituting the 
integrands from Eqs. ®-([TT]) furnishes 

f" 1 _ ^ 

(l+p2)3/2 - - — 

+ Nrr - Ntt - P{rp + u)w' = 0, 

{Nrrw'rpY -I- P{rp -I- m)(1 -I- u') = 0 
with boundary conditions 

^glRo^^^^d — ^rri^RoO^RoO ^ 

7gii?o sm0d = Nrr{Roo)w{Roo)Roo at z = h, 

Isg - Isl = Igl cos Oc at .S = 0 . 

Here, dc and 6d are the contact and dihedral angles that the interface Fg; forms 
with the plate and membrane, respectively. The last boundary condition only 
holds when the meniscus is a bridge, while the boundary term associated with 
r]i(£) when the meniscus is a bubble is always zero with r'{£) —> oo. Equa¬ 
tion (ITTl) is the Young-Laplace equation, and Eqs. (1^ and (1^ are the in- 
and out-of-plane equations of equilibrium for membranes in von Karman’s the¬ 
ory with the capillary pressure acting normal to the neutral plane. Note that 
Eqs. (13011 and m demand w'{Roo)P{h) = 1, implying w' ^ oo as 9d ^ 7 r/ 2 , 
contradicting the assumptions of von Karman’s theory [d^. To resolve this 
issue, we undertake a scaling analysis in section [3] to simplify Eqs. (l30l) and 

(ED). 


(27) 

(28) 

(29) 

(30) 

(31) 

(32) 
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Figure 2: Schematic of perturbations to (a) simply connected (bubble), and (b) 
doubly connected (bridge) menisci. 


2.3 Stability from second variation 

The Jacobian matrices of the functional integrants in Eqs. (|6]) and © are not 
symmetric. Since analyzing the second variation for functionals of multiple func¬ 
tions with non-symmetric Jacobians is intractable 0 , we simplify the problem 
by neglecting in-plane variations at the hole edge, prescribing (j)i{Roo) = 0. 
Accordingly, the second variation is 


2d^Et ^ 

£ J Zo ^ 2 

dz+ f (drp 

^ Roo 


JRoo ^ 


(33) 

where ^ = F — PK, ^ = G — PM, and 


c^if) _ ® 

fp'j'l ^ 

dz 

(34) 

agi-a.) —ca,, 

£)(u) _c^ _, 

■=^ — ^UU t ^uu' 5 

dv p 

(35) 

agiw) 

^ — ’^W'W' 1 

aiw) _ 

- ^WW 1 ^WW' ^ 

dvp 

(36) 

= 

^ TV ' 

r 

(37) 

Note that the boundary term at zq 

in Eq. (1331) arises only when the meniscus 


is a bubble, so S^Et has only one boundary term at h when the meniscus is 
a bridge. Furthermore, when the meniscus is a bubble, is not bounded due 
to the axial symmetry condition at Zq = i (see Appendix]^, and it is unsuit¬ 
able for representing quadratic forms. To resolve this issue, using the mapping 
(z,^i) —>■ {y,N), the first integral and boundary terms in Eq. (l33l) are repre¬ 
sented with respect to the normal variation N{y), where y = cos6 (see Fig.j^la) 
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and Appendix . Equation (1551) then reduces to 


Q = Igi f [(1 - + ‘^yNN' - N^] dy + - igiN^l) 

J I/n 2/0 


+ 


'Vo 

f-fl r 


' Roo 


c 


rpC(j)f + ( — + Pw' ) (^1 


^ drp + / {Nrr + Cw'‘^)rp'ijj'^drp, (38) 
-Roo 


+ 


nh _ 


Q = 


dz + 


^0 

' Roo 


rpCcj)'^ +(—+Pw' 


rR 

^Roo 


drp + / {Nrr + Cw'‘^)rp'tp'^drp, (39) 


for simply connected and doubly connected menisci, respectively. Here, Q = 
25'^Etle'^, 2/0 = cos 9q corresponding to the polar angle at the hole edge (Fig.[5Da)), 
and 


= — 


IgP 


IgV 


p/(l _j_ p/2^3/2 


(1 +p/2)3/2’ 


^{r) ^ _ 


Igl 


r(l + 


(40) 


A necessary (sufficient) condition for Et to have a minimum is that Q be 
non-negative (strongly positive) for an equilibrium solution Il8| . Here, strong 
positivity, referred to as nonlinear stability in the literature l45ll , must be dis¬ 
tinguished from positive-definiteness. According to Vogel [^, Q > 0 does 
not imply a strict local minimum for constrained infinite-dimensional problems. 
Nevertheless, defining stable equilibria as those for which Q > 0 holds, Mad- 
docks [ 2 ^ derived sufficient conditions for the positive-definiteness of the second 
variation. Interestingly, Vogel ^ showed that these conditions are sufficient for 
Madoccks’ functional to have a strict minimum. Therefore, following Maddocks 
[23 |. we adopt Madoccks’ definition of stability to relate stability to the slope 
of equilibrium branches. 

Our analysis differs from Madoccks’ theory in two respects: (i) The elasto- 
capillary energy Et is a functional of multiple functions, the stationary points 
of which are represented by non-smooth functions. Moreover, Et has three sep¬ 
arate spectra, which complicates the relationship between stability exchanges 
and equilibrium-branch turning points, (ii) The boundary condition at the hole 
edge, where the membrane and meniscus meet, is not prescribed. Here, pertur¬ 
bations are arbitrarily finite, posing the same difficulties as discussed by Myshkis 
et al. [ 2 ^ and Vogel for analyzing the stability of capillary surfaces with 
free contact lines. 

The quadratic forms in Eqs. (|58ll and (l39|) demand ^ ^ 

where perturbations satisfy the boundary conditions Eqs. (|T6|-([2n|). Since each 
bilinear term in Eqs. (1551) and ((551) is bounded, Q can be represented as 0,0 


Q = Igl {■^(f)N,N)^ + 


(41) 
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Q = (.i?(f)Ci,Ci)i + {^(u)(t>i,4>i)-^ + (■^(^)'0i,V'i)i (42) 

for simply connected and doubly connected menisci, where (•, •)j^ is the inner 
product [l| and , ^[f) j -^{u) j -^{w) are uniquely determined by the respective 
bilinear terms in Eqs. (1381) and (1391) . Because perturbations are arbitrary at 
the hole edge, the spaces from which perturbations are selected are not 
symmetric, and, thus, are not generally self-adjoint. On 

the other hand, integrating Eqs. (I55|) and (15^ by parts leads to 

Q = jgi + (.if(*)'0i,'0i) , (43) 

Q = (44) 

for simply connected and doubly connected menisci, where -S?(y) = U"--^ — PJ”y) 
with double primes denoting the second Frechet derivative and perturbations 
satisfying the boundary conditions 


N'{vo) - N{yo)/yo = 0, 7V(1) = finite, 

$((o) = 0, eUM + = 0, 

<Pi{Poo) = 0, 4>i{R) = Oj 

'0'i(-Roo) = 0, '(/’i(^)=0, 

furnishing a symmetric for evaluating Q. Here, 


^ N = - — 
dy 


dy 


- 2iV, 


(45) 

(46) 

(47) 

(48) 


(49) 


= - 


“ dr„ 


dz 


dz 






d(t)i 
p J 


C 


+ ( — + Pw' ) «f>i, 


~ dr. 


[Kr+Cw''^)'^ 

ar„ 


'p 


(50) 

(51) 

(52) 


subject to Eqs. (H5]) - (H51) are regular Sturm-Liouville operators; thus, they are 
self-adjoint and Fredholm Establishing a relationship between the spec¬ 

trum of the barred operators on and those of unbarred operators on 
significantly simplihes the analysis, providing a setting to apply the well-studied 
Sturm-Liouville theory. Lemma 2.5 of Vogel [i^ furnishes this relationship by 
stating that the barred operators have the same number of negative and non¬ 
positive eigenvalues as the corresponding unbarred operators. Moreover, the 
constant-volume constraint {vi = const.) can be written 


since 5vi = 5J. 
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3 Scaling analysis 


Starting from Eq. (|28)) and taking Rqq <C R, we begin by scaling all lengths with 
R. Given NrrjC = u' + w''^f2 + vuffp and Ntt/C = uu' + l2-\-ulrp, where 
u = u/R^ w = w/R, and fp = Vp/R, we find u ^ viP' io balance all the terms in 
the in-plane equilibrium. Considering the bending and stretc hing parts of the 
strain energy in von Karman’s theory for axisymmetric plates |40l | 




w 


/2 



2vw'w'‘ 


Tpdrp, 


Q.S 




u w 


2vuu' 

Tp 


lyuw'^ 

Tp 



Tpdrp, 


and noting that Vp ^ R, w ^ H, C ^ Eb, and D ~ Eb^, where D and b are 
the bending rigidity and plate thickness, one infers fls ^ when b/H 
1 by comparing the energy scales Hb ~ DE[^/R^ and fls ^ CH'^/R^. For 
thin membranes, this justifies neglecting the bending energy compared to the 
stretching energy. 

Given Nrr ~ C{w/R)^, |w| ~ {\Qc\/ follows from the out-of-plane 
equilibrium, where k = Rqo/R, Qc = PRm/lgii Nc = Cjlgi are the scaled 
hole radius, scaled capillary pressure and elastocapillary number. Noting that 
|?i;| -I-kA = n, the elastocapillary number corresponding to a specific state of the 
system in Fig. [T]can be estimated as Nc ~ \Qc\/— kA)^, where A = h/Roo 
and n = El/R are the meniscus slenderness and aspect ratio. For example, 
when the meniscus is a bubble contacting the plate, before bridging the gap, as 
will be discussed elsewhere, Qc = —2sm9d and A = cot(0c;/2); thus, at 9^ = 9c, 
which is the critical dihedral angle below which the bubble cannot be stably 
pinned to the hole edge [2^, we have 


Nc 


2 sin 9c 

K(n — Kcot 9c/2)^ 


(54) 


From Eq. (l30l) . Nrr{Roo) ~ Igi, while Nrr ~ Cw'^ in the main body of the 
membrane. Thus, Nrr{Roo) ^ Nrr when Nc ^ implying that Eqs. (150]) 

and (1311) can be approximated by a free-edge boundary condition for slender 
cavities when Nc is large. 


4 Membrane profile 

The stretching part of the strain energy leads to nonlinear equations of equilib¬ 
rium, which can be solved numerically in most practical problems. Although 
the case considered in this paper, namely, axisymmetric plate with a hole at the 
center, has a series solution |40l |. expressions for the unknown coefficients are 
cumbersome. Therefore, we apply a variational approximation, as commonly 
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(a) (b) 



Figure 3: Comparison of the numerically exact solution (solid) and variational 
approximation (dashed) of (a) in-plane and (b) out-of-plane displacements: 
Qc = —2, K = 0.1, v = 0.3, and Nc = 100,500,1000, 5000,10000 (downward in 
(a) and upward in (b)). 



r^/R t^/R 

Figure 4: Comparison of the numerically exact solution (solid) and variational 
approximation (dashed) of axial forces in the (a) radial and (b) tangential di¬ 
rections: Qc = —2, K = 0.1, ly = 0.3, and Nc = 100,500,1000,5000,10000 
(downward). 
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used in the literature [6|,l25|, to construct a general solution for the membrane 
deflection. 

Here, we approximate the boundary condition at Vp = i?oo as a free edge. 
Noting that the test function 


w = woCl - ^p) 


(55) 


is consistent with the boundary conditions and the plate stress distribution 
under tension at zero deflection zho = {Qc/NcY^^Ku, is obtained by mini¬ 
mizing the stretching energy, where 





i[l - V + K^(l + v)] 
k ‘^)(1 — — ly + k 2(1 + ly )] 


(56) 


which is in agreement with the scaling relation for w derived in section [3] Equa¬ 
tion ([53 contrasts with the test function w = zZ’o(l ~ fpY that Mastrangelo & 
Hsu [23 used to describe the bending and stretching contribution to the overall 
deflection of beams. The approximation accuracy hinges on choosing appropri¬ 
ate test functions for bending- and stretching-dominated regimes. Figures|3|and 
[4] demonstrate a reasonable agreement between the variational approximation 
and numerically exact solutions of the membrane equilibrium, for a wide range 
of deflections. 


5 Second variation spectra 

We adopt the foregoing variational approximation in this section to determine 
the second variation spectra corresponding to the in-plane and out-of-plane 
equilibria. 


5.1 In-plane spectrum 

The spectrum of the membrane in-plane equilibrium is determined by 


r ^(^u)Z = rpdZ 
\ ZiRoo) = 0, Z(i?) = 0, 


where z?, Z, and rp are the eigenvalue, eigenfunction, and weight function of 
■ Non-dimensionalizing Eq. (EH) furnishes 


"-i; YY:) + - ‘I = “ 

Z{k) = 0, Z{1) = 0, 


(58) 


with i?* = and B = 2{Q(.lNcY/^K^I k. The general solution of Eq. (157)) 

is Z{fp) = C'iJi(mfp) -I- C 2 Yi(mfp) with m = \/B where Ji and Yi are 
the Bessel functions of first and second kind. A similar eigenvalue problem was 
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Figure 5: In-plane spectrum i?* with i = 0,1, 2, 3 (upward). 




derived by Timoshenko & Gere for the buckling of circular plates under 
in-plane compressive loads. Imposing the boundary conditions furnishes 




Jl(m^rp) - - - rYi{mirp) 

Yi{miK) 


i = 0,1,2,• 

where mi is the ith root of 


i(TO)Yi(mK) — Ji(TOK)Yi(m) = 0. 


(59) 

(60) 

(61) 


From Eq. dSQl), unless B = ttiq at a given k, Zq corresponding to dp = 0 
has only a trivial solution where ker(.5f({i)) = {0} and zi = 0 at P = 0. Here, 
0 is the identically zero function, and the overdot denotes differentiation along 
equilibrium branches. Therefore, stability loss due to in-plane perturbations is 
not generally related to pressure turning points. Studying these instabilities, 
which are responsible for wrinkling in thin elastic membranes IJ, Il5l . l3l| . is 
beyond the scope of the present work and will not be further elaborated upon. 

Figure [5] shows the first four eigenvalues of the in-plane spectrum. Note 
that mg Ri 16 is accurate in the range 0 < k < 0.2 to within 12 % of computed 
values. A positive spectrum is ensured by i? < mg. As water is removed from 
the elastocapillary system of Fig. [U B reaches its maximum value when the 
membrane touches the plate where Bmax = /[K^k{1 — k^)^]. To ensure 

that "d* > 0 always holds during drying, we require Bmax < ^g, leading to 


n < 23/4(1 _ ^2), 


(62) 
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based on the foregoing approximation of mg and the scaling relation ^ I/k. 
Consequently, Nc ^ guarantees that the in-plane spectrum 

is positive, and that the membrane profile can be accurately approximated by 

Eq. (1551) . 


5.2 Out-of-plane spectrum 

The spectrum of the membrane out-of-plane equilibrium is given by 

r if(^)r = rpAT 
\ Y'{Roo) = 0, ZiR) = 0, 


(63) 


where A, Y, and Vp are the eigenvalue, eigenfunction, and weight function of 
Without attempting to solve Eq. (l63)) . we demonstrate that has a 
positive spectrum. From the quadratic form 

r,) = A, {rpY,,Y^ = / rp{N,r + Cw'^)Y^^drp, 

J Rqq 



it follows that Ai > 0 because Nrr > 0, based on the variational approximation 
discussed in section 0] Similarly to the in-plane spectrum, Yq corresponding to 
Ao = 0 has only a trivial solution where ker(.if(^)) = {0} and w = 0 aX P = 0. 


5.3 Meniscus spectrum 

When the meniscus is a bridge, the spectrum cannot be determined analytically. 
Therefore, in this section, we only study the meniscus spectrum for simply 
connected menisci determined by 

/ ( 04 ') 

1 yo^'(yo) -= 0, X(l)= finite, ^ 

where fj, and X are the eigenvalue!! and eigenfunction of ). These also denote 
the eigenvalue and eigenfunction of for doubly-connect menisci. Solving 
Eq. (IMl) furnishes 

W.(y) = Ci,,P„,(2/), (65) 

/r* = [(2m. + 1)2 - 9]/4, i = 0,l,2,--- (66) 

where Pm is the real-valued order Legendre function of first kind, and rm is the 
fth root of 

yoPm. (yo) - Pm. (yo) = 0. (67) 

The boundary condition of Eq. (I64|) at y = yo, resulting from perturbations 
that can displace the hole edge, is a key feature of this study. It implies that the 
meniscus contact line at the hole edge exhibits a mixed characteristics of free and 
pinned contact lines, depending on the bubble size. In the limit yo —>■ 0, where 

®Not to be confused with the chemical potential in Eqs. Q and Q. 
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(a) (6) 




Figure 6: Meniscus spectrum; (a) order of eigenfunctions rrii and (b) corre¬ 
sponding eigenvalues with i = 0,1, 2,3 (upward). 


the bubble is a hemisphere, the contact line behaves similarly to a pinned contact 
line. As shown in Fig. [51 this limit, which corresponds to the pressure turning 
point (P = 0) of the elastocapillary model in Fig. [1] with simply connected 
menisci, occurs at the point of stability exchange where /xq = 0. Therefore, 
as expected, this limit coincides with the stability limit of pressure-controlled 
spherical menisci with a pinned contact line 27|. Here, the order of the Legendre 
function takes integer values where rui = 2z -|- 1, and ker(Jf(i!)) ^ {0}. An 
important implication is that, unlike the in-plane and out-of-plane spectrum, 
f 0 at P = 0. 


6 Stability along equilibrium branches 


In this section, a relation between stability and the slope of equilibrium branches 
in vi versus P diagrams for constrained and unconstrained problems is estab¬ 
lished. We assume that the foregoing variational approximation for the mem¬ 
brane equilibrium and, particularly, Eq. (1621) always hold, so > 0 and > 0 
for all i. As discussed in section [2731 stability can be determined by the sign 
of Q in Eqs. (H51) and with N, satisfying Eqs. (H5l) - (H51) . 

Following the Ritz method 18| , Q is examined in a countable dense subspace of 
the existence of which is guaranteed by the separability of Q, because 
there always exits a function in the dense subspace that is arbitrarily close to 
any / £ Furthermore, the eigenfunctions of the Sturm-Liouville opera¬ 
tors in Eqs. (|^ - (I5^ form a complete orthogonal basis in and, thus, 

span the respective dense subspace. Hence, (j)i G span{Zi}, ipi £ spanjFi}, 
N £ spanjAi} when the meniscus is a bubble, and £ span{Xi} when the 
meniscus is a bridge. 

We unify the representation of quadratic forms in this section by expressing 
Q in terms of for simply connected and doubly connected menisci. When the 
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meniscus is a bubble, because (^(f)CijCi) = Igi {■^{r)N, N), the sign of Q is 
determined by the eigenvalues of .5f(f), no matter how Q is expressed. 

Lemma 1. The necessary condition for an equilibrium branch to be stable in 
the unconstrained problem is fj-i > 0,'di > 0, Ai > 0 for all i. 

We assume the contrary is true. For example, let /tq < 0 and pi > 0 for i > 1. 
Therefore, po < 0 . Choosing = oq^o with 


ao > 




1/2 


IMoI (-’^ 0 , ^q)i 


results in Q < 0, which is a contradiction. Here, p and X denote eigenvalues 
and eigenfunctions of ). □ 

Note that in lemma 1, Eqs. (|4T]) - (|42]) are used to illustrate how the relation¬ 
ship between the spectrum of barred and unbarred operators can be applied 
to prove stability. Hereafter, Eqs. (I43ll - (l44l) are directly used to determine the 
sign of Q. The following lemma connects stability to the slope of equilibrium 
branches in the unconstrained problem. 

Lemma 2. The slope of a stable equilibrium branch at any point in the J versus 
P diagram is non-negative. 

Differentiating Eqs. (l24)l - (|M)l along an equilibrium branch results in 

= P'Kv) ’ y = r,u,w, ( 68 ) 

furnishing 


P'i'r) - PJ(f)]fp) + {{U'i'u) - PJiu)K^) + - PJiu,)H 


= P 


J[f) X) + { J(u) I ■“ / + ( '^(w) I 


'(u) 




Using Eq. 


j =U[f),r 


fKr' -\- zK 


1 h 


+ ( J[u) + 


+ ( J[w) I ) + 


uMu 


WMy, 


R 


Roo 
1 R 


Rlh 


-\- RoRoh -l- 
J Roo 2 


Substituting for K and M from Eqs. (0 and dlB, all the boundary terms cancel 
each other, leading to 

j = 


and, consequently. 


([[/(")-P 4 )]A^) + ([c/(C) 





= pj = p^ 


dP' 
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Since the branch is stable, fXi > 0, > 0, Ai > 0 according to Lemma 1. Thus, 

the left-hand side is non-negative. □ 

Similarly, the following lemma connects stability to the slope of equilibrium 
branches in the constrained problem. 

Lemma 3. Suppose that, on a segment of an equilibrium branch, is non¬ 
singular and has precisely one negative eigenvalue. Then, the segment is stable 
if and only if the slope at any point in the J versus P diagram is negative. 

We prove the lemma for doubly connected menisci. Consider the following 
perturbation decompositions 

h='Vr + ar]r, 4'i = Vu + ariu, 'ipi = + arj^,, 


where 

= J[y)^ y = r,u,w. (69) 

Because are all non-singular, ker(.5f(j;)) = {0} and € ker(.if(y))-*-. 

Therefore, rjy always have a solution because are Fredholm operators. From 
the volume constraint in Eq. (1531) . 


T 'Cu) T {^•^{w)Vw: ‘Cw') 


furnishing 


(X Vu) “b j Vw')'\ • 

Note that (.if(ii)77„, ? 7 „) , ?7i«) , > 0 because 

the in-plane and out-of-plane spectrum are positive. We first show that the 
necessary condition for Q > 0 is 

T Pr') “b “b (,<^(^w)'ejw j Pw') ^ O- 

We assume the contrary holds. Choosing Vr = oqXo with 

r/ \ / \ 5 

'Cu/ F 

|^oKNo,Afo) 

results in Q < 0, which is a contradiction. Therefore, it is always possible to 
construct perturbations that lead to instability if Eq. (1701) does not hold. Here, 
^0 and Xq are the negative eigenvalue and corresponding eigenfunction of . 
Next, we show that Eq. dZO]) is sufficient for Q > 0. Since perturbations are 
selected from a countable dense space, rjr can be written rjr = + biXi 

such that 

“I 1/2 

, i^it ^’^{u)Vu-> Vu') “t“ ■) Vw) 

ImoI (Afo, Afo) 
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(a) (6) 



Figure 7: Equilibrium branch of the elastocapillary model in Fig. [T]for simply 
connected menisci with k = 0.1,1/ = 0.3, Nc = 15000,11 = 0.2; (a) numerical 
computation and (b) schematic representation of pressure and volume turning 
points. 


for Eq. ((701) to hold. Moreover, any arbitrary meniscus perturbation can be 
written = aoXg + Choosing a = ao/bg leads to Vr = 

implying that I'r) > 0, and, consequently, Q > 0. Hence, all meniscus 

perturbations, including those with (^i,Xo) ^ 0, can be decomposed into rjr 
and Vr, leading to a strictly positive second variation, provided Eq. (701) holds. 
Substituting Eq. (|69|) into Eq. (l68)) furnishes 

Vy = y/P^ y = r,u,w. 


giving 


Vr') + {■^(u)Vui Vu) + (,-^(w)Vw t Vw) 
1 ■ 

~ P - 


J[f) I ^ ) + ( 






j _ dJ 
P~dP’ 


which completes the proof. □ 

The proof of lemma 3 for simply connected menisci similarly proceeds using 
the decomposition N = Vm + ctVrn and accounting for the relation between 
N and given by Eq. (A.4() . Figure [7] shows how lemmas 1-3 can be applied 
to determine stability from the shape of an equilibrium branch in the uncon¬ 
strained and constrained problems. In the unconstrained problem, where P is 
the control parameter, only the segments AB and CE can be stable according 
to lemma 2. The stability of the segment AB, excluding B, is deduced from the 
system configuration at A, corresponding to the fully saturated state. Here, the 
meniscus and membrane are planar, and the system is evidently stable, imply¬ 
ing that iJi, Ai > 0 for all i. Assuming that varies continuously along the 
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(a) (&) (c) 



Figure 8 : Meniscus equilibrium (solid) and perturbed (dashed) states at vol¬ 
ume and pressure turning points in Fig. [T] The most dangerous perturba¬ 
tion normalized by {N,N) = Rqq at (a) the pressure turning point B where 
N{y) = •\/3i?ooPi(2/) and (b) the volume turning point C where N{y) = 
y/T^^/ (Pmo^Pmo) RooPmoiy)- (c) A Safe perturbation N{y) = CoXo{y) + 
CiXi{y) at C normalized by {N,N) + ('(/'i,'0i) = Rqq where 4‘i{rp) = 0 and 
tjji{rp) = ao-l-aifp-|-a 2 rp such that Eqs. (fT6)) - (l20|) and (|5^ are satisfied, leading 
to Q > 0 . 


equilibrium branch, stability is lost at i?, and the entire segment BCDEf is 
unstable. The segment AB is also stable in the constrained problem, where vi 
is the control parameter. Moreover, the stability of BC, along which has 
one negative eigenvalue, is deduced from lemma 3. Beyond the volume turning 
point C, the entire segment CDEf is unstable with respect to constant-volume 
perturbations. 

Determining the stability of equilibrium branches is essential for predicting 
the dry-state conformation of the present model. This is illustrated by an ex¬ 
ample in Fig. [3 Here, the point D corresponds to the state where the bubble 
is tangent to the plate. Further decrease in vi forces the bubble to bridge the 
membrane and plate, which is a necessary step for the elastocapillary system of 
Fig. [T]to collapse. However, when drying from the fully saturated state at A, 
collapse does not occur for the given parameters in Fig. (3 because the system 
loses stability at C before the bubble can bridge the membrane and plate. 

Note that spherical menisci with a pinned contact line are always stable to 
constant-volume perturbations [1^. Furthermore, as discussed in section [3 the 
membrane is stable for all deflections, provided Eq. (1621) is satisfied. Therefore, 
the meniscus and membrane are individually stable along the entire branch 
ABODEf. However, the elastocapillary system as a whole subject to vi = 
const, is unstable along CDEf, revealing an intimate connection between stabil¬ 
ity and the coupling of elastic and capillary forces. This manifests in the bound¬ 
ary shared by the elastic and capillary part where the meniscus and membrane 
interact through boundary displacing perturbations. Moreover, the nature of 
instabilities are influenced by the control parameter and how the meniscus and 
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membrane interact with each other, as demonstrated in Fig. [S] 


7 Concluding remarks 


We have developed an elastocapillary model to study drying-induced structural 
failures, such as those arising from stiction in microelectromechanical systems. 
The model comprises an elastic membrane and a meniscus, deformed by the 
same pressure differential, interacting through a shared boundary. The exis¬ 
tence of a stable equilibrium branch from the fully saturated to collapsed state 
is an essential precursor for structural failures. We examined the model stability 
and equilibrium using variational and spectral methods. Stability was related 
to the slope of equilibrium branches in the liquid content versus pressure dia¬ 
gram for the constrained and unconstrained problems. A variational approxi¬ 
mation, complemented by scaling analysis, was derived, furnishing closed-form 
expressions for membrane equilibria. This approximation leads to a positive 
out-of-plane spectrum. For a given geometry, there is a critical elstocapillary 
number above (below) which the in-plane spectrum is positive (has a negative 
eigenvalue). These in-plane instabilities are a common cause of wrinkling in thin 
membranes. Thus, except for thin membranes, only meniscus perturbations can 
be dangerous for the elstocapillary system. This paper extends the work of Mad- 
docks [ 2 ^ to elastocapillary systems that are subjected to boundary displacing 
perturbations, revealing a close connection between stability and the coupling 
of elastic and capillary forces. We demonstrated that the stability of the menis¬ 
cus and membrane alone does not imply that the elastocapillary system as a 
whole is stable; the destabilizing effect of the membrane and meniscus interact¬ 
ing through their shared boundary must also be accounted for. Moreover, our 
results support a general concept in catastrophe theory that stability exchanges 
occur at turning points with respect to the control parameter, thereby, reducing 
costly stability computations to search for folds on equilibrium branches. 


Appendix A 


Expressing perturbations in is problematic for axisymmetric simply connected 
menisci because ^ L^. Assuming that perturbed states are also axisymmetric 
(j.e., dz/dr = 0 at z = £), we have 


dz _ 1 -I- 77 (£ H- 

dr f' -I- -!-••• 


1 + 771 f' 


?/2 


■e -I- = 0 at £ = 


(A.1) 


in view of Eqs. (1^ and (|T^ . Since f'{i) oo and {£) = finite, it 

follows that 1^1 (•^) I ^ O(f') —>■ oo satisfies the axial symmetry condition. More¬ 
over, Erom Eq. (1^ . |^i(£)| —oo. Therefore, f is unbounded, and, thus, unsuit¬ 
able for representing perturbations when the meniscus is a bubble. Here, rep¬ 
resenting functional variations with respect to the normal variations of menisci 
resolves the issue. 
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When the meniscus is a bubble (sphere), it is convenient to express variables 
as functions of the meridian-curve arclength s or the polar angle 9 = s/Rg (see 
Fig. [2]). The displacement vector from the meniscus equilibrium states to its 
perturbed states is written 

dx , ,. , 

— =^ir-f?7iz. (A.2) 

£ 

Given h = sindf — cos0z, the normal variation 


(5x 

iV = h • — = sin 0 — ? 7 i cos 9 
e 


(A.3) 


is obtained, furnishing 


6 


N 

sind’ 


(A.4) 


resulting in dA^/ds = —dyyi/ds at z = i. Moreover, rj'i —>• ^^/r' and d? 7 i/ds —>■ 
as z ^ i, giving d? 7 i/ds(^) — >■ 0 even if |,^(| ~ 0{f'). Therefore, 
dN/ds = 0 and N = finite at z = I, implying N £ L^. 

The following formulas are useful for representing the meridian curve when 
the meniscus is a bubble: 


furnishing 


f = \/Rj- {zc - z)^ = Rs sin 9, 

r' = = cotd, 

r 

f" = _i±i^ = 1 

f Rs sin^ 9 ’ 


(A.5) 

(A.6) 

(A.7) 


= Rsjgi sin^ 9, 



COS 6^ 


(A.8) 


where Zc and Rs are the z-coordinate at the center and radius of the sphere. 


Appendix B 


Consider the functional 

Ay\= / y: [xa,xb] r, (b.i) 

of continuously differentiable functions y defined on a variable domain where 
the branches of stationary points are parametrized with t, and the stationary 
points are represented by y = y{x,t). Then, differentiating the functional and 
its first Frechet derivative along a branch furnishes 


j = 



yFy> -I- xF 


(B.2) 
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(B.3) 


where 



dx 
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